Probing the Priedmann equation during recombination with future CMB experiments 
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We show that by combining measurements of the temperature and polarization anisotropies of the 
Cosmic Microwave Background (CMB), future experiments will tightly constrain the expansion rate 
of the universe during recombination. A change in the expansion rate modifies the way in which the 
recombination of hydrogen proceeds, altering the shape of the acoustic peaks and the level of CMB 
polarization. The proposed test is similar in spirit to the examination of abundances of light elements 
produced during Big Bang Nucleosynthesis and it constitutes a way to study possible departures 
from standard recombination. For simplicity we parametrize the change in the Friedmann equation 
by changing the gravitational constant G. The main effect on the temperature power spectrum is 
a change in the degree of damping of the acoustic peaks on small angular scales. The effect can 
be compensated by a change in the shape of the primordial power spectrum. We show that this 
degeneracy between the expansion rate and the primordial spectrum can be broken by measuring 
CMB polarization. In particular we show that the MAP satellite could obtain a constraint for the 
expansion rate H during recombination of SH/H ~ 0.09 or SG/G ~ 0.18 after observing for four 
years, whereas Planck could obtain SH/H < 0.014 or SG/G < 0.028 within two years, even after 
allowing for further freedom in the shape of the power spectrum of primordial fluctuations. 



I. INTRODUCTION 

The parameters of our cosmological model will be de- 
termined with great accuracy by upcoming data from 
Cosmic Microwave Background experiments, galaxy sur- 
veys, weak lensing surveys, Lyman alpha forest studies, 
and other observations. With the new data it will be 
possible to perform a number of consistency checks that 
will strengthen our confidence in the underlying model. 
Some of these consistency checks have already been per- 
formed with existing data. Recent analysis of the CMB 
data have resulted in constraints on the baryon density 
Slbh 2 that are in excellent agreement with its determi- 
nation based on the study of the primordial abundances 
of light elements (e.g. The combination of CMB 

data with local measures of the Hubble constant 0] and 
measures of the local strength of galaxy clustering result 
in a determination of the cosmological constant that is in 
good agreement with results from the study of the lumi- 
nosity of distant supernovae (e.g. Recently also 
joint BBN and CMB constraints on different dark energy 
models have been discussed in [(J. 

One of the assumptions of the cosmological model that 
has been hard to test is the explicit validity of the Fried- 
mann equation - the relation between the expansion rate 
of the universe and its matter content. The difficulty 
lies in finding an epoch in the evolution of the universe 
during which both the energy density and the expansion 
rate can be determined independently. The two obvious 
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candidates are the present time and Big Bang Nucleosyn- 
thesis (BBN). 

Precise measures of both the expansion rate and the 
matter density in the local universe arc difficult to accom- 
plish and suffer from various systematic problems. At 
present the expansion rate has been measured with er- 
rors on the order of 10% 0. Direct determinations of the 
matter density however are more uncertain. It is fair to 
say that the general conclusion from these studies is that 
the Friedmann equation only holds true if either a cos- 
mological constant or a curvature term is added. This is 
because direct determinations of the present matter den- 
sity almost always point to Q m < 1 (e.g. see |f|). Neither 
the cosmological constant nor the curvature scale can be 
constrained independently, that is without going through 
the Friedmann equation, so at best we can say that the 
Friedmann equation has not been tested accurately at the 
present epoch. A more radical interpretation would be 
that the Friedmann equation has been tested but that 
the test has failed. Although we do not support this 
interpretation, at present some infrared modification of 
gravity cannot be ruled out observationally. Such modi- 
fications are being explored for example as ways to solve 
the cosmological constant problem (see for example Q) 
or to explain the accelerated expansion inferred from the 
luminosity distance to high redshift supernovae without 
resorting in a cosmological constant or a quintessence 
field (e.g. see 00). 

During the epoch of Big Bang Nucleosynthesis (BBN) 
the situation is more fortunate. The energy density is 
dominated by radiation which we think we can estimate 
accurately. On the other hand the expansion rate affects 
the freeze-out abundances of light elements, so that a pre- 
cision test of the Friedmann equation can be performed. 
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The standard procedure is to constrain the number of rel- 
ativistic degrees of freedom, , which can be translated 
into a limit on the number of neutrino species. A lot of 
progress has been made in determining the primordial 
abundances. Recently the deuterium abundance in hy- 
drogen clouds at high redshift was accurately determined 
[ToL 111! Il2| . Building on a prior of N„ > 3 these data 
have been exploited to enforce an upper limit of 3.2 at 
2cT for the number of neutrino species, based purely on 
BBN considerations [l3|. 

Progress has also been made in appreciating the sys- 
tematical uncertainties that impair the determination of 
primordial 4 He (see e.g. Built on the safely es- 

tablished abundance ranges for Deuterium, Helium and 
Lithium, it can be shown that the uncertainty in the 
number of relativistic degrees of freedom during BBN is 
around 9 % (68% C.L.)[l||. This constraint can equally 
well be phrased as a constraint on the validity of the 
Friedmann equation: during Nucleosynthesis the ratio of 
the squared expansion rate to the energy density can de- 
part by only 9% from what is predicted by the Friedmann 
equation. Constraints on the expansion history during 
BBN have also been established in [l7|. 

In this paper we propose using the anisotropics in the 
CMB to perform a test similar to the one that has been 
done using BBN. We will show that such a test can ulti- 
mately constrain the validity of the Friedmann equation 
during recombination more accurately than what has so 
far been reached within Nucleosynthesis, albeit in a more 
model dependent way. 

During recombination, the energy density is dominated 
by the density of non-relativistic matter which we can- 
not estimate directly. However the dark matter energy 
density enters in two different ways and one can exploit 
this to simultaneously determine the dark matter density 
and the expansion rate during recombination. 

The ratio of matter to radiation energy density sets 
the redshift of matter radiation equality. At that time 
the expansion rate changes from a scaling as t x l 2 to t 2 / 3 . 
Perturbation modes of the photon-baryon fluid that en- 
tered the horizon during the radiation dominated era be- 
have differently than those that entered during the mat- 
ter dominated era. Modes that entered during radia- 
tion domination provided the dominant contribution to 
the total density perturbation that generated the gravita- 
tional potential. On the contrary modes that entered the 
horizon during matter domination, were sub-dominant in 
their contribution to the total density perturbation which 
was dominated by the dark matter fluctuations. The 
gravitational potential acts as a source for perturbations 
in the photon baryon fluid. As a result, small scale modes 
that entered the horizon in the radiation era go through 
a sort of feedback loop that increases their amplitude as 
they cross the horizon (for a review of CMB physics see 
for example 18j). The anisotropy power spectrum is very 
sensitive to the redshift of matter radiation equality and 
thus the CMB should very accurately determine the ra- 
tio of dark matter to radiation energy density, i.e. the 



parameter Q m h 2 . 

The dark matter density dominates over other energy 
components in the Friedmann equation during recombi- 
nation. In the standard scenario, it sets both the red- 
shift at which matter and radiation become equal and 
the rate of expansion during recombination through the 
Friedmann equation. In this paper we break the link be- 
tween energy density and expansion rate by introducing 
a free parameter to modify the Friedmann equation. We 
investigate how well such a parameter can be constrained. 

From a pragmatic perspective our study can be re- 
garded as the investigation of a particular departure from 
standard recombination. Different variations have been 
studied in the literature. For instance the possibility that 
energetic sources of Ly-a-photons could be present dur- 
ing recombination to delay it was considered in |l9j. Also, 
the possiblility of a time variation of the fine structure 
constant was investigated [U |H In [23| the ef- 

fects of a time dependence of the gravitational constant 
have been outlined. The conclusion of these investiga- 
tions and ours is that with future CMB data departures 
from standard recombination will be severely constrained 
and perhaps modifications that point to interesting new 
physics could be discovered. 

In section^ we will introduce our model and in section 
lllll we will investigate the constraints that can be set with 
currently available data and forecast what future CMB 
experiments might be able to achieve. We will conclude 
in section llVl with discussion. 



II. THE MODEL: VARIATION OF THE 
GRAVITATIONAL CONSTANT 

In this section we introduce the model we will use to 
investigate how well one can test the Friedmann equation 
using the CMB. The problem is somewhat more subtle 
than in the case of BBN because we are dealing with 
the dynamics of perturbations that could be affected by 
the "new physics" in ways other than through a change in 
the expansion rate. Thus we need to find a self-consistent 
way of modifying both the dynamics of the universe and 
that of the perturbations. 

One possibility is to add another component that con- 
tributes to the energy density during recombination, in- 
creasing the rate of expansion at that time. The addi- 
tional component could be a quintessence field with a 
potential and initial conditions tuned so that it has some 
effect during recombination and is unimportant or only 
marginally important at other times (except perhaps to- 
day when it could start to dominate) . This approach has 
the virtue of only modifying the expansion rate but it 
introduces too much freedom because results depend on 
when exactly this extra component is important. In such 
a model we also expect the ratio of the sound horizon 
at recombination to the angular diameter distance to the 
last scattering surface to change and thus that the acous- 
tic peaks be slightly shifted. At late times the evolution 
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of the gravitational potentials will also induce an inte- 
grated Sachs- Wolfe (ISW) effect. As a result, constraints 
on any specific model of this kind will come from these 
three effects [24|. 

In our study we want to isolate the information en- 
coded in the change of the expansion rate at recombi- 
nation so we will use a simpler prescription and assume 
that the gravitational constant G is somewhat different 
from its locally measured value. We introduce a single 
parameter A such that, 



G -> A G. 



(1) 



The expansion rate is proportional to A. With this pre- 
scription not only the Friedmann equation gets modified 
but also the dynamics of the perturbations changes be- 
cause it depends on the strength of gravity. We will show 
in the following that our prescription has the nice feature 
that it only changes the CMB power spectrum through 
the change in recombination, allowing us to isolate the 
observable effects of this change. The basic reason is 
that gravity does not have a preferred scale and that 
we only measure angles when studying the CMB. If G 
were slightly different all that would happen is that the 
universe would be expanding a bit faster or slower by a 
factor A so that the "expansion clock" would be running 
at a different rate. Such a change cancels in the ratios of 
distances that we measure with the CMB. The only way 
we can find out that such an alteration had occurred is 
by having an independent clock that measures the ex- 
pansion rate. In our case this independent clock will be 
the physics of hydrogen recombination. In this sense our 
simple test is very similar to what has been done in the 
context of BBN. 



A. Effect of A on the CMB anisotropies 

The dependence of the Hubble parameter and the dy- 
namics of perturbations on the gravitational constant will 
lead to modifications of the CMB anisotropies as we vary 
the parameter A. We will discuss the physics in this sec- 
tion. 

We start by considering the modification to the Fried- 
mann equation, 



H 2 



8tt 



Gp 



8tt 



VGp 



(2) 



where p is the total energy density. As a function of the 
expansion factor a and A, the expansion rate H satisfies: 



H(a,X)=Xf(a) 



(3) 



where the function /(a) is independent of A. Thus with 
this simple prescription, the shape of the function H of 
a is not changed by A, only the amplitude changes. For 
example the redshift at which matter and radiation con- 
tribute equally to the energy density docs not change. 



The change introduced is a simple rescaling of the "ex- 
pansion rate clock" . 

In order to understand how the anisotropies get mod- 
ified, we start by writing down the integral solution for 
the temperature anisotropies produced by a mode of 
wavevector k observed towards direction n |25| . The tem- 
perature can be written as an integral along the line of 
sight over sources, 



AT(n,k) = / dr S(k,T)e i: 



k-nD(r) 



9(r) (4) 



In this equation S(k, r) is the source term, g(j) is the 
visibility function, and D(t) is the distance from the ob- 
server to a point along the line of sight corresponding to 
the conformal time t (adr = dt). Dots indicate differen- 
tiation with respect to t. 

The visibility function g(r) can be written in terms of 
the opacity for Thomson scattering k as 



g(r) = Kexp(— n) = —d/drexp^—n) 



with 



K = (7t 



an e (r)dr, 



(5) 



(G) 



where ut is the Thomson scattering cross section and 
n e (r) is the number density of free electrons. We have 
also defined k = aTan e . Finally, the source term in the 
integral equation is given by 



n I °7 

S = cp+ — + n ■ v b 



(7) 



where <fi is the gravitational potential, 5 7 is the fractional 
perturbation in the photon energy density and Vf, is the 
baryon velocity. 

The acoustic oscillations in the photon-baryon plasma 
satisfy (see e.g. |2g) 



&y + 



R 



{1 + R) 



R 



-6- -k 2 (j) 
R v 3 v 



(8) 



with the sound speed c| = 1/3(1 + R), and the baryon- 
photon momentum density ratio R = (pt, + pf,)/(p^ + 
p 7 ) ~ 3/9fc/4p 7 . The velocity satisfies, 



5 1 + kv.y + = 



(9) 



Finally the gravitational potential satisfies the Poisson 
equation 



-k z <f> = 4ir\ 2 Gp5 total , 



(10) 



where p$ total gives the combined perturbation due to all 
the fluids. 
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We are now ready to study the dependence of AT on 
A. For this purpose it is best to consider the expansion 
factor as a time variable rather than r. We note that, 

d da d d , , d , 

— = —— = a ■ H ■ — = A/(a)a — . 11 
ar ar da aa aa 

As a result, when we change time variables, every time 
derivative introduces a factor of A. By inspection of equa- 
tions (JSJl and llOjl it is clear that the dynamics of a mode 
with wavcnumber A: in a universe with A ^ 1 is equivalent 
to the dynamics of a mode with k! = k/X in a universe 
with A = 1. That is, 



S(k,a,X) = S(k/X,a,X = 1). 



(12) 



We have explicitly included the A dependence of the 
source to make our argument clearer. 

To obtain the CMB power spectrum, we first need to 
expand equation i n Legendre polynomials. The am- 
plitude of the I expansion coefficient is 



AT;(fc, A) 



da S(k, a, X)ji{kD{a, A))<?(a, A). (13) 



We have introduced g(a, A) = —d/daexp(—K,). The con- 
formal distance D is given by 



D(a, A) 



da 



a H{a) a i 



X~ 1 D(a, X = 1). 



(14) 



Thus if the visibility function where to be independent 
of A we would have, 

ATj(fc, A) = ATi(fc/A, A = 1). (15) 

The power spectrum is calculated from ATJ(fc, A) using, 

CiiX) = J T P(k)\A Tl (k,X)\ 2 



dk 
~k 



* 7 P(k'X)\A Tl (k',X = l)\ 2 , (16) 



where P(k) is the power spectrum of primordial fluctu- 
ations which is usually taken to be a power law P(k) cx 
k n ~ x . Thus we see that provided we adjust the am- 
plitude of the primordial power spectrum appropriately 
Q(A) = Q(A = 1). 

Our result is qualitatively very easy to understand: 
gravity introduces no preferred scale, so the dynamics 
of the perturbations remains the same when scales are 
measured in units of the expansion time. As a result, the 
angular power spectrum does not change as we change A. 

Of course this conclusion only holds true if the visibility 
function is not affected by A. However the physics of re- 
combination does introduce a preferred timescale, so the 
power spectra of the anisotropics will actually change. In 
other words, in our simple minded prescription the only 
source of change is the difference in the way recombi- 
nation proceeds as we change the expansion rate of the 
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FIG. 1: Ionization fraction as a function of redshift for three 
values of A = 0.25, 1,4. 



universe at recombination. This is the sense in which our 
model resembles the studies done in the context of Big 
Bang Nucleosynthesis. 

Let us now turn to study how the visibility function 
changes with A. It depends on the ionization fraction 
x e = n e /n,H, where n e again is the free electron density 
and iih is the number density of hydrogen atoms. The 
evolution of the ionization fraction is modified when G is 
changed. It evolves according to (e.g. pT^l: 



dx e 

— — = aC r 
dr 



!3{T b ){l-x e )~n H a^(T b )xl] (17) 



where a(t) is just the scale factor, f3{Tb) is the collisional 
ionization rate from the ground state and (T&) is the 
recombination rate to excited states. The baryon temper- 
ature is Tt and the Peebles correction coefficient (which 
also depends on the expansion rate) is denoted C r [28| . 
The transformation from r to a(r) as a time variable 
using equation 1|11[1 . makes clear, that contrary to what 
happens to the perturbation equations, x e (a) depends on 
A. We plotted x e for different values of A in Figure^ The 
behavior is easy to understand; the faster the universe is 
expanding at a given redshift (i.e. the larger the A), the 
more difficult it is for hydrogen to recombine and hence 
the larger is x e . 

The change in x e leads to a change in the visibility 
function which we show in Figure [21 As A is increased, 
the visibility function becomes broader. This broadening 
leads to a larger damping of the anisotropies on small 
(angular) scales, as shown in Figure [31 We note however 
that even for a factor of four change in A the changes in 
the visibility function are rather small. What happens 
is that if a given redshift after the 

start of recombination increases. However, when calcu- 
lating optical depths this change almost exactly cancels 
with the decrease in the time intervals between different 
redshifts due to the increased expansion rate. As a result 
changes in both the location and shape of the visibility 
function are small even for large changes in A. 
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FIG. 2: The visibility function as a function of conformal time 
for A = 0.25, 1, 4. The axes have been rescaled to take out the 
overall scaling of r with A. 



Figure |21 shows that the effect of A is to change the rel- 
ative amplitudes of the acoustic peaks on different scales. 
This effect can be compensated by changing the relative 
amplitude of modes of different scales in the primordial 
power spectrum. 

We are now going to study what happens to CMB po- 
larization and to show that it can lift the degeneracy be- 
tween A and the shape of the primordial power spectrum. 
To understand the effect of A on the polarization we will 
employ the simple analytic expression for the amplitude 
of the Q Stokes parameter produced by a single Fourier- 
mode k El: 



Q oc CskSro sm(fccgTD)e 



k /k D 



(18) 



where td is the conformal time corresponding to the peak 
of the visibility function, Stjj is its width and hp de- 
scribes the damping of the small scale modes. As we dis- 
cussed above, the damping increases with the width of 
the visibility function, so the exponential factor will lead 
to the same effect we described for the temperature. The 
difference in the case of the polarization is the extra Std 
in the amplitude of the polarization. This extra factor 
comes from the fact that if the visibility function is wider 
the photons will travel on average longer between their 
last two scatterings which will enhance the quadrupole 
anisotropy and will thus lead to a higher level of polar- 
ization [23 . For that reason, there exists a characteristic 
wavemode value k*, for scales larger than which the po- 
larization power spectrum will increase with A, while it 
will decrease for scales smaller than k* . 

The polarization power spectrum is plotted in Figure 
0J We see that for Is larger than around 800, the polariza- 
tion behaves just as the temperature does, decreasing for 
increasing A. On larger scales the effect is opposite, the 
amplitude of polarization relative to temperature roughly 
increases by 10% when A increases by 20%. This response 
of the anisotropics on A is what will help to break the de- 
generacy between A and the primordial power spectrum 



when information from polarization is included. 

For completeness we show the temperature polariza- 
tion cross correlation power spectrum in Figure [SJ The 
cross correlation will be easier to detect in experiments 
such as the MAP satellite where the accuracy of the po- 
larization measurement is limited by detector noise. The 
cross correlation power spectrum behaves similarly to the 
polarization power spectrum; when A is increased the 
power on small scales is suppressed while on large scales 
it is amplified. 



III. CONSTRAINTS ON A 

In our likelihood analysis of currently available and 
simulated future data we let vary Wdark matter = 

^dark matter h 2 , Wbaryon = ^baryori^ 2 , the Optical depth 

due to reionization r r i, A and the amplitude and shape 
of the primordial power spectrum. We explicitly assume 
that the universe is flat. 

As we mention above, we expect there to be a degener- 
acy between the shape of the primordial power spectrum 
and the parameter A. To study this effect in the case of 
future satellite missions which will measure polarization 
and could break this degeneracy, we introduce additional 
freedom in the shape of the spectrum. Rather than just 
assuming that P(k) is of power law form 



P(k) = k 11 - 1 & 



lnP(fe) 
ln(fc) 



= 71—1, 



(19) 



we also allowed spectra with curvature by adding another 
term in the expansion of In P as a function of Ink. 

\nP(k)/P(k ) = (n-l)\n(k/k )+a[\n(k/k )} 2 + ... (20) 

where fco is the pivot point. With this prescription the ef- 
fective slope of the power spectrum changes slightly with 
scale, 



d1nP(k) 
d\nk 



= n — 1 + a\n{k/ko). 



(21) 



We employed two kinds of likelihood analysis. For an 
evaluation of what currently available data can tell us 
about A we used an importance sampling Markov-chain 
method to generate a large number of cosmological mod- 
els distributed according to the likelihood distribution 
£(model|data). This is more efficient than an exploration 
of the entire parameter space, because the sampling is 
weighted and statistics can be established over the target 
distribution itself. Moreover the algorithm is very easy to 
implement. The weighted sampling is achieved through 
the Metropolis-Hastings algorithm. The application of 
the Markov method to the extraction of cosmological pa- 
rameters from CMB information has been suggested in 

To speed up the power spectrum computations for the 
Markov chain we made use of the k-splitting technique 
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FIG. 3: Effect of A on the Temperature power spectrum: Higher peaks are more severely damped as A increases while the 
height of the first peak is almost unchanged. 




discussed in |31| . For each model we marginalized ana- 
lytically over the amplitude of the scalar fluctuations. In 
the end we constructed histograms and computed expec- 
tation values and variances for each parameter directly 
from the Markov chain. 

In order to estimate what the satellite missions MAP 
and Planck will be able to tell us about the relation be- 
tween the energy density and the expansion rate of the 
universe during recombination, we investigated the shape 
of the likelihood function C(Ci\0i) (for a power spectrum 
Ci given a model consisting of the cosmological param- 
eters 9i) in the vicinity of its maximum directly by a 
Fisher matrix evaluation. This method has been widely 
used to make predictions for the errors that are to be ex- 



pected in the extraction of cosmological parameters from 
planned CMB-experiments. The Fisher matrix is given 
by the expectation value of the second derivative of the 
logarithm of the likelihood function C(Ci\9i). Assuming 
Gaussianity of the likelihood it is of the form 

l A.B ° V i 

where A and B run over the three observables: tem- 
perature, E-type polarization, T-E cross correlation and 
i,j run over the cosmological parameters. The covari- 
ance matrix between parameters is given by the inverse 
of the Fisher matrix. Overall we verified a good agree- 
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FIG. 5: The effect of A on the cross correlation between temperature and E-type polarization. 



Par am. 


mean 


""Markov 


A 


1.749 


0.471 


ns 


1.038 


0.0553 



TABLE I: Mean value and standard deviation a for the cur- 
rently available temperature data. 



merit between Fisher matrix and Markov chain results 
which confirms that the likelihood function C(Ci\9i) re- 
sembles relatively well a Gaussian in the vicinity of its 
maximum value. 



A. Constraints from currently available 
temperature data 

In order to find the constraints which can be imposed 
on the parameter A using available temperature data, we 
employed a compilation of 30 experiments which func- 
tions as a complete account of pre-MAP CMB tempera- 
ture anisotropy information |42|, The data, which have 
been compressed to 29 bins, are plotted in Figure |SJ In 
this context, we actually assumed a simple power law 
behaviour of the primordial scalar perturbations. 

The results we obtained from the Markov-chain anal- 
ysis are shown in Table [I] and Figure As expected 
from the rather weak dependence of the anisotropies on 
A we find that current data cannot put severe constraints 
on the expansion rate during recombination even though 
other cosmological parameters are well determined. 

Having emphasized the importance of measuring the 
linearly polarized component of the CMB, we should also 
note that its recent first detection by DASI [3(| , unfortu- 
nately has too large errorbars to deliver much informa- 
tion about A. Adding the DASI data in fact only changes 
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FIG. 6: Compilation of current CMB temperature anisotropy 
data, superimposed are three models with differing values of 
A. 




FIG. 7: Histogram for likelihood of A given present CMB 
temperature anisotropy information. 
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the standard error in the determination of A by 3%. 



As a fiducial model we have adopted the parameter 
values 



B. Future satellite missions and the relevance of 
measuring polarization 

An up to date estimate of the expected angular reso- 
lutions and sensitivities of the MAP and Planck satel- 
lites has been obtained from the experimental groups 
websites. For both satellites we combined the three fre- 
quency channels with the highest angular resolution and 
took into account the number of polarized instruments. 
In the case of MAP the sensitivity estimate has been 
provided for 2 years of observation and the angular reso- 
lutions (9f w hm are 13.2, 21.0 and 31.8 arcminutes for the 
three channels observing at 90, 60 and 40 GHz. This 
leads to a raw sensitivity of 

Wp 1 = (O.IU^K) 2 

If MAP observes for four years the raw sensitivities (w" 1 ) 
will be halved. For the Planck satellite mission (here the 
sensitivity estimates are given for a one year observation 
period), the three channels (217, 143 and 100 GHz) at 
#fwhm = 5.0,7.1 and 9.2 arcminutes give together 



w,, 



(0.0084/^ 
(0.0200/xFsT) 5 



Again the raw sensitivities for a two year observation will 
be half of these values. For both satellites, a sky coverage 
of fsky = 0.8 was assumed. 

From these experimental characterisics the full estima- 
tor covariancc matrices for each multipole I can be con- 
structed (e.g. |37j). The diagonal terms of the covariance 
matrices for temperature, polarization and cross correla- 
tion are 



Cov(C Tl ) = 



(Cti + w T 1 B l 2 ) 2 



(21 + l)f S ky 

Cov(C 2 El ) = 2 (Cei + w E l B t 



Cav(CS,) 



(21 + l)f sky 
1 



(21 + l)f Sk y 

x (C E i + w P 1 B 



(23) 
(24) 

[Oct + (Cti + w T x B{ 2 ) (25) 



(26) 



where the beam window function Bi is to be constructed 
from the relevant frequency channels with their individ- 
ual sensitivities w r as 



Bf = 



-l{l+l)8 2 b ,c 



(27) 
(28) 



Here, the standard width of the beam 9b is obtained from 
the full width half maximum resolution by 



■'fwhn 



V8~h72 



(29) 



r 


Q K 


n A 






n 3 


A 


0.05 





0.7 


0.14 


0.02 


1 


1 



(30) 



They imply a hubble constant of h = 0.73. 



1. Expected constraints from MAP 

The analysis of expected constraints from the MAP 
satellite mission shows the constraints it can put on A are 
rather weak. Better sensitivity is needed to accurately 
determine the polarization and higher angular resolution 
to map the damping tail. On the other hand our re- 
sults illustrate how the inclusion of a curvature term for 
the primordial spectrum significantly weakens the con- 
straints that one obtains from the temperature data. 

A precision test of the Friedmann equation will have 
to wait until both polarization and the damping tail are 
measured accurately. In the near future experiments such 
as Boomerang, BICEP, Polatron and others are expected 
to significantly improve polarization measurements while 
CBI, ACBAR and others will map the damping tail. 
In a few more years, the Planck satellite will measure 
both temperature and polarization accurately enough to 
severely constrain any change in A. We will present an 
analysis of Planck's sensitivity in the next section. 

In Section [H] we discussed the different response to 
a change of A of the temperature and the polarization 
anisotropy. We found that on large angular scales po- 
larization power is increased if we increase the expan- 
sion rate during recombination, while the temperature 
anisotropy is almost not affected on these scales. When 
only temperature is being measured, changes in the ex- 
pansion rate of the universe during recombination are 
strongly degenerate with the slope of the spectrum of the 
initial scalar perturbations, which can be clearly seen in 
Figures OH and [5] Because polarization responds differ- 
ently on different scales to a change of A, it breaks this 
degeneracy. We show this in Figure The result of the 
likelihood analysis which includes polarization informa- 
tion is included in Tabic ITT1 

We noticed that a major part of the information on A 
from MAP's polarization measurements will come from 
the cross correlation between temperature and polariza- 
tion. When adding just the cross correlation information 
we found that we gain almost 90 % of the information 
that is gained in the case in which all three estimators 
are included. 



2. Expected constraints from Planck 

The Planck satellite explores the very small structures 
in the primeaval plasma in a multipole range up to nearly 
/ = 3000. Even if only Planck's temperature data are 
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FIG. 8: Likelihood distributions of A for a one year observation with MAP. Compared is the case in which one assumes a simple 
power law spectrum of the initial perturbations (solid lines) with the case in which one leaves further freedom to P(k) (dashed 
lines). To obtain the Figure on the left, only temperature information was used and for the one on the right polarization was 
included. 
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FIG. 9: We plot the 68.3, 95.4 and 99.7% contour lines of 10 5 cosmological models of our Markov chains. The upper two plots 
show a projection onto the A — a-plane; the left plots are for temperature and the right ones include polarization information. 
The lower two plots show the degeneracy of A with the slope of the spectrum of scalar perturbations, ns- It can be seen in the 
upper and lower plots that polarization information reduces the degeneracy with the shape of the initial power spectrum. 
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T 




(F^ /2 , a 


T+P 




(F^f /2 , a 


A 


0.1928 (0.1279) 


0.2616 (0.1992) 


A 


0.1373 (0.0903) 


0.1658 (0.1165) 


n 


0.0294 (0.0248) 


0.0646 (0.0595) 


n 


0.0170 (0.0144) 


0.0489 (0.0418) 


a 


X 


0.0237 (0.0222) 


a 


X 


0.0151 (0.0137) 



TABLE II: Fisher matrix results for MAP's expected 2 year data (in brackets are the values expected for a 4 year observation). 
The first two columns are the results when polarization information is not included. In the second column the curvature of the 
primordial spectrum was left to vary. In the last two columns polarization information was included. 
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(F^) 1/2 


(F^) 1/2 , a 


T+P 




(F^ 2 , a 


A 


0.0170 (0.0152) 


0.0325 (0.0278) 
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0.0115 (0.0093) 


0.0174 (0.0141) 
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0.0118 (0.0106) 


0.0182 (0.0160) 
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0.0072 (0.0060) 


0.0098 (0.0080) 
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0.0072 (0.0620) 
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0.0039 (0.0033) 



TABLE III: Fisher matrix results for Planck's estimated 1 year data (brackets contain the 2 year results). Columns are 
equivalent to those in Table ITT1 
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used, the standard error in A will after a one year ob- 
servation be as small as 3.2%, even if one marginalizes 
over the parameters that describe the shape of the pri- 
mordial power spectrum fTable lTTTl) . This corresponds to 
a constraint on the gravitational constant of = 0.064. 
Finally, the improvement gained from Planck's polariza- 
tion data ~ 0.017 <^> ~ 0.034) is shown on the 
right hand side of that same Table. If one assumes that 
there is no curvature in the primordial power spectrum 
the constraint on G would be as small as ^ = 1.8% after 
an observation of two years. FigureHUlandllllshow again, 
analogously to the case of MAP, how polarization helps 
break the degeneracy between A and the parameters that 
describe the shape of primodial perturbations. 

We have extended our investigation beyond Planck to 
the case of an experiment which has essentially no noise 
and where the errors in the CMB power spectra are on all 
scales dominated by the cosmic variance term. For exam- 
ple experiments currently being considered for measuring 
the B-modes of CMB polarization would be cosmic vari- 
ance limited for E polarization over a wide range of Is 
(e.g. Such an optimal experiment, exploring struc- 

tures into a multipole range of I — 4000, represents the 
limit of how much information on A one could extract 
from the CMB in principle. We found an expected error 
for A of order 0.3% which translates into a constraint of 
the value of the gravitational constant during recombi- 
nation of ~ 0.6%. 

IV. DISCUSSION 

We have made the expansion rate of the universe a 
free parameter in a likelihood analysis within an eight- 
dimensional cosmological parameter space. For simplic- 
ity we assumed that the gravitational constant is changed 
by a factor A. We showed that an increase of A leads 
to a wider visibility function which in turn increases the 
damping of anisotropies on small scales and increases the 
level of large scale polarization. 

We calculated the constraints that current CMB data 
can impose on the expansion rate. The constraints that 
can be imposed on the parameter A using the information 
from the damping tail are severely weakened by our lack 
of knowledge about the shape of the primordial power 
spectrum. We showed that measuring polarization helps 
to break this degeneracy. Current data can only con- 
strain A to about 47% at 1 a. We showed that MAP 
could obtain 9% error bars for A while for Planck error- 
bars go down to 0.9%. We also explored the ultimate 



limit that could be achieved by a cosmic variance limited 
experiment measuring anisotropies up to I — 4000 and 
found errors of under a percent in that case. Thus next 
generation experiments should be able to deliver very ac- 
curate constraints on the expansion rate of the universe 
during recombination. 

We acknowledge that if the variation of the gravita- 
tional constant during recombination is taken seriously 
a model needs to be built where G changes after re- 
combination and converges towards the stable value ob- 
served in laboratory experiments today and where its cur- 
rent rate of change is less than the experimental bound 
G _ io- 12 yr -i[23|, [39|. If we introduce a scalar field 
to control the value of G we would also have to require 
that this field does not lead to an unacceptably large fifth 
force and does not violate solar system constraints such 
as shiftin g th e orbit of the moon through the Nordvedt 
effect Hfllll. 

The shift of G after recombination will induce a change 
in the angular diameter distance to recombination, shift- 
ing the CMB power spectrum in I. We have shown that 
future experiments will be able to constrain the change of 
G to a few percent due to its effect at recombination. As 
a result the induced shift in the peak positions would be 
small and could be interpreted as slightly different val- 
ues of Q m and/or Q\, the two parameters that control 
the distance to the last scattering surface in conventional 
models. In the same way, any induced integrated Sachs- 
Wolfe (ISW) effect would be difficult to observe because 
of the cosmic variance limitation. If observed with a high 
enough accuracy it should not be identical with what is 
predicted by a simple cosmological constant model. 

In this paper we have shown that future measurements 
of the CMB anisotropy will be able to extract informa- 
tion about the relation between the expansion rate and 
the energy density of the universe during recombination, 
because of its effect on the recombination history of hy- 
drogen. 
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